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The purpose of this study was to evaluate the effect of rectangular region of 
interest (ROJ) size on modulation transfer function (MTF), to develop the 
radial ROI, and to compare both ROIs performances for MTF measurement 
using a cylindrical polymethyl methacrylate (PMMA) phantom. The PMMA 
phantom used in this study was rotated 45°. Four rectangular ROIs and a radial 


ROI were created to measure the MTF value. The rectangular ROI sizes were 


3x41, 21x41, 41x41, and 61x41 pixels; each was placed at upper phantom 
edge. The radial ROI’s length was 41 pixels and placed at several points in 
phantom edge. The MTF calculation was automatically conducted using 
MATLAB. The MTFs from rectangular ROIs and radial ROI were then 
compared. The comparison of the MTF measurement was also conducted 
using three different filters. The MTF which used radial ROI was smoother 
Modulation transfer function than those of rectangular ROI for all filters. This indicated that radial ROI was 
Polymethy! methacrylate more resistant to noise than rectangular ROI. Rectangular ROI with the 41x41 
phantom pixels had similar 50% and 10% MTF values with the radial ROI. The MTF 
value which was obtained using radial ROI is more accurate and robust than 
those obtained using rectangular ROI. 
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1. INTRODUCTION 

The quality of computed tomography (CT) image needs to be maintained by conducting regular 
quality control (QC). Spatial resolution is one of the important parameters measured in QC, which was 
conducted in a metric of modulation transfer function (MTF). The MTF calculation is initially introduced by 
Judy [1], with the main purpose was to identify frequency domain response of imaging systems [2], [3]. MTF 
measurement was also used to characterize small structure [4], [5] and to compare the performance of imaging 
systems [6]. Thus, accurate MTF measurement is needed. 

MTF can be measured using point, line and edge objects where the objects are spread and lead to the 
point spread function (PSF), line spread function (LSF), and edge spread function (ESF) [7], [8]. The MTF 
would then be obtained by applying Fourier transform to the LSF and PSF which would produce | -dimensional 
(1D) and 2-dimensional (2D) MTF, respectively. Furthermore, MTF could also obtained by differentiating and 
applying Fourier transformation to the ESF [9]-[11]. 

Kayugawa et al. [12] conducted a study to find out an optimal size for square region of interest (ROT) 
used for calculating accurate and precise MTF using CatPhan 600 and a thin wire phantom. The results showed 
that the optimal sizes for square ROI were 50 pixels for FC10 filter, 38 pixels for FC50 filter, and 40 pixels for 
FCS52 filter when using bead method. Optimal sizes for square ROI were 42 pixels for FC10 filter, 39 pixels 
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for FCS50 filter and 44 pixels for FC52 filter when using wire method. Therefore, the size of the ROI used in 
the process of MTF calculation can affect the accuracy of the MTF. There are several phantoms used to 
practically measure MTF, such as CatPhan 600 [12], [13], ACR [14], [15], wire [16], [17] and AAPM CT 
performance phantom [18], [19]. In developing countries, those phantoms were only owned by few hospitals. 
Recently, there were studies that conducted MTF measurement using PMMA phantom [20]. Anam et al. [20] 
proposed a MTF measurement using the PMMA phantom. The PMMA phantom was rotated 45 degrees 
(clockwise/counter-clockwise) to avoid the ionization pencil holes. This MTF measurement was compared 
with using wire phantoms. The difference of MTF measurement using the upper of edge PMMA and wire 
phantoms was only within + 4% [20]. However, in that study, a square ROI at circular edge of phantom was 
used and it may decrease the MTF value because circular phantom edge is linearly averaged in x-axis direction. 
Therefore, in the current study, the effect of rectangular ROI size on MTF measurement will be analyzed. We 
hypothesize that decreasing ROI size will yield a more accurate MTF. However, small ROI size may be 
susceptible to the noise. In addition, we will develop radial ROI which will be placed at several point in edge 
of PMMA phantom and compare the resulted MTF to those from rectangular ROIs. 


2. METHOD 
2.1. ROI determination 

The head PMMA phantom with a diameter of 16 cm was used in this study. The phantom was scanned 
with a Toshiba Alexion-4 CT scanner with three different filters: FC13, FC30, and FC52. The parameters of 
scanning were field of view (FOV) of 20 cm, 120 kVp, 100 mA, rotation time of | s, slice thickness of 2 mm, 
and using axial mode. In this study, the phantom was rotated 45 degrees. In the phantom, there are five 
ionization chamber pencil holes. Figure 1 shows the position of the holes in each image, with filter FC13 
(Figure 1(a)), FC30 (Figure 1(b)), and FC52 (Figure 1(c)). 

In this study, the algorithm from Droege and Morin [19] was improved using MATLAB software. 
Previously, Anam et al. [20] used 41x41 pixels ROI and the measurement was performed automatically. 
Figure 2 shows various ROI size used in this study. The current study used four sizes of rectangular ROIs: 3x41 
(Figure 2(a)), 21x41 (Figure 2(b)), 41x41 (Figure 2(c)), and 61x41 pixels (Figure 2(d)). The radial ROI was 
automatically measured along the edge of phantom (Figure 2(e)). Each line of radial ROI had length of 41 pixels. 

The radial ROI was formed with four steps. First, the phantom image was segmented. Second, the 
center of the image was determined. Third, a line which passed through the edge of phantom image 
perpendicularly was yielded. The last, ROIs one by one along the edge of phantom image in 35 to 55 degrees 
except in the ionization pencil holes were formed. 


2.2. MTF measurement 

MTF measurement was performed automatically, which was started by taking all profiles of the ROIs. 
For rectangular ROIs, the value of each pixel along the x-axis was averaged to yield the ESF. For radial ROIs, 
the profile of each ROI was taken and then averaged to yield single averaged profile of ESF. Each ESF was 
interpolated to place four additional data. The interpolated ESFs were then differentiated, zeroed, and 
normalized to yield the line spread function (LSF) curve. Fourier transform was then applied to the LSF curve 
in order to obtain the MTF curve. The value of MTF at 10% and 50% were automatically determined. The 
steps of MTF determination using radial ROI could be seen in the Figure 3. 


(b) 


Figure 1. Images of phantom with different filters were rotated 45 degrees; Filter (a) FC13, (b) FC30, and 
(c) FC52 
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(d) (e) 


Figure 2. The different ROIs were used for measuring MTF on head PMMA phantom. Four rectangular ROI 
sizes are (a) 3x41 pixels, (b) 21x41 pixels, (c) 41x41 pixels, (d) 61x41 pixels, and (e) radial ROI 
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Figure 3. The steps of MTF determination using radial ROI 


3. RESULTS AND DISCUSSION 

Figure 4 shows the results of image filtered with FC13 from all ROIs, with it’s ESF (Figure 4(a)), 
LSF (Figure 4(b)), and MTF (Figure 4(c)). The curves of the ESF is shown to be overlapping with one another, 
indicating that all MTFs are comparable. However, the ESF curve of ROI size of 3x41 pixels is shown to be 
more fluctuating due to the image noise, as predicted. These lead to a more fluctuating LSF curve because of 
the differentiation process. As a result, the MTF curve produced from this ROI also became more fluctuated. 
On the other hand, the measurement which used ROI size of 61x41 pixels produce the lowest peak at LSF 
curve. The MTF curve produced by this ROI was also decreasing faster than the other, indicating that it 
produces lower spatial resolution, as predicted. The 10% and 50% MTF values are tabulated in the Table 1. 

Figure 5 shows the results of image filtered with FC30 from all ROIs, with it’s ESF (Figure 5(a)), 
LSF (Figure 5(b)), and MTF (Figure 5(c)). It can be seen that all ESF curves for radial and all rectangular ROI 
sizes have similar form, except for 3x41 ROI. Both tails of the LSF curves for all ROIs were also shown to 
have negative peaks. In the MTF curve, radial ROI has the highest peak of + 1.2. The MTF curve of 3x41 
rectangular ROI obtained by using this filter was also shown to be fluctuating. The 10% and 50% MTF values 
are tabulated in the Table 2. 
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Figure 4. The comparison of (a) ESF, (b) LSF, and (c) MTF using FC13 filter for rectangular ROI sizes of 
3x41, 21x41, 41x41, 61x41 pixels, and radial ROI 


Table 1. The values of 10% and 50% MTFs for the rectangular ROIs and radial ROIs using FC13 filter 


Spatial resolution (cycles/mm) 


ROI types ROI size (pixels) 


MTF 10% MTF 50% 
3x41 0.71+0.05 0.38+0.04 
Reotanzular 21x41 0.71+0.01 0.38+40.01 
= 41x41 0.69+0.01 0.3740.00 
61x41 0.62+0.01 0.3340.00 
Radial 1x41 0.69+0.00 0.37+0.00 
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Figure 5. The comparison of (a) ESF, (b) LSF, and (c) MTF using FC30 filter for rectangular ROI sizes of 
3x41, 21x41, 41x41, 61x41 pixels, and radial ROI 
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Table 2. The values of 10% and 50% MTFs for the rectangular ROIs and radial ROIs using FC30 filter 
ROI types ROI size (pixels) —_ oe pine 


3x41 1.11+0.07 0.66+0.22 

Rectangular 21x41 1.07+0.02 0.7240.02 
41x41 1.01+0.03 0.70+0.01 

61x41 0.86+0.03 0.60+0.01 

Radial 1x41 0.93+0.00 0.70+0.01 


Figure 6 shows the results of image filtered with FC52 from all ROIs, with it’s ESF (Figure 6(a)), 
LSF (Figure 6(b)), and MTF (Figure 6(c)). In the ESF, LSF, and MTF curves, the rectangular ROI size of the 
3x41 pixels was also shown to have large fluctuation. The MTF of 61x41 pixels produces lowest spatial 
resolution. The 10% and 50% MTF are tabulated in the Table 3. 
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Figure 6. The comparison of (a) ESF, (b) LSF, and (c) MTF using FC52 filter for rectangular ROI size of 
3x41, 21x41, 41x41, 61x41 pixels, and radial ROI 


Table 3. The values of 10% and 50% MTFs for the rectangular ROIs and radial ROIs using FCS2 filter 
ROI types ROI size (pixels) eee gas a 


3x41 1.55+0.64 0.8740.12 

Rectansular 21x41 1.11+0.03 0.81+0.03 
41x41 1.07+0.02 0.77+40.01 

61x41 0.94+0.02 0.65+0.01 

Radial 1x41 0.95+0.01 0.76+0.00 
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Figure 7 shows the MTF comparison for 41x41 pixels and radial ROI for all filters, with FC13 shown 
in Figure 7(a), FC30 in Figure 7(b), and FC52 in Figure 7(c). The curves which was using FC13 filter were 
shown to coincide. The value of 10% and 50% MTF for radial ROI are slightly higher than 41x41 pixels of 
rectangular ROI. However, MTF curves using FC30 filter for radial ROI is smoother than 41x41 pixels of 
rectangular ROI, where 50% MTF of radial ROI is slightly higher than 41x41 rectangular ROI but 10% MTF 
of radial ROI is lower than 41 x 41 rectangular ROI. MTF curves using FC52 filter for 41x41 rectangular ROI 
is slightly higher than radial ROI, but the overall curve of radial ROI is slightly smoother than 41x41 
rectangular ROI. The 10% and 50% MTFs as shown in Table 4. 
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Figure 7. The MTF comparisons for rectangular ROI of 41x41 pixels and radial ROI using (a) FC13, 
(b) FC30, and (c) FC52 filters 


Table 4. The values of 10% and 50% MTFs for the 41x41 pixels of rectangular ROI and radial ROIs using 
FC13, FC30, and FCS2 filters 
Spatial resolution (cycles/mm) 
ROI type FC13 FC30 FC52 
MTF 10% ~MTF50% MTF10%  MTF50% MTF10% MTF50% 
Rectangular (41x41) 0.6940.01 0.37+40.00 =1.01+0.03 = 0.7040.01 =: 1.07+0.02 _—_(0.77+0.01 
Radial 0.69+0.00 _0.3740.00 _0.93+0.00 _0.70+0.01 —_0.95+0.01 _0.76+0.00 


This study developed the previous algorithm for automated MTF calculation on head PMMA phantom 
[20]. In the current study, 3 filters were investigated and produce different spatial resolution. The size and 
shape of the ROI were also shown to affect the accuracy of the MTF. Radial ROI has many lines of ROI and 
they can be used as additional representation of ESF and to decrease error due to noise [21]-[23]. The ESFs of 
the radial ROIs are averaged to produce single ESF [24]. 

The MTF curves from rectangular ROIs and radial ROI using FC13 filter were coinciding with each 
other, indicating that the MTFs are comparable. The rectangular ROI sizes of 3x41 pixels and 61x41 pixels for 
all filters are slightly different. The ROI size of 3x41 pixels has greater fluctuation than those with wider ROIs. 
The decrease of the ROI size resulted in less number of data points used to calculate the MTF. The smaller 
number of data is more affected by image noise [25]. Larger noise caused the MTF to have more fluctuations [22]. 
Fewer data can produce bias error caused by noise. Whereas, the rectangular ROI size of 61x41 pixels leads to 
lower spatial resolution for all filters due to the curved edge within the ROI is linearly averaged. This indicates 
that wider ROI will lead to lower spatial resolution. 10% and 50% MTF values using rectangular ROI sizes of 
41x41 pixels, however, were similar to the radial MTF values for all filters. The radial ROI was created from 
several points in phantom edge [24]. The radial ROI for various filters was more resistant to noise due to many 
lines which were used for MTF calculation. Therefore, it can be used to accurately measure MTF on PMMA 
phantom which has circular edge. 

The radial ROI method is accurate as long as the spatial resolution of the CT does not pass the Nyquist 
frequency. The Nyquist frequency of this imaging condition is 1.28 mm’. It can be seen from Figure 7(c), that 
the MTF which was using radial ROI dropped faster than those which used rectangular ROI because the 
frequency was close to the Nyquist frequency. If the spatial resolution of the CT system is greater that Nyquist 
frequency, smaller FOV should be used. However, that condition will cause the image to be truncated. 
Therefore, the MTF measurement using radial ROI requires alternative approach, such as using PMMA 
phantom with diameter of 8 cm and/or scanned with FOV of 12 cm or less. The 16 cm head PMMA phantom 
scanned with 24 cm FOV can also yield oversampling MTF using rectangular ROI method, if shifting and 
rebinning data approach are applied [25]. These approaches will be investigated in the future study. 
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4. CONCLUSION 

The comparison of radial and rectangular ROIs has been successfully conducted. The results of this 
study show that the size of rectangular ROI affects accuracy of the MTF. It is found that MTF values of 
rectangular ROI with size of 41x41 pixels close to those from the radial ROI. Therefore, the optimal rectangular 
ROI size for MTF measurement using edge of PMMA phantom is 41x41 pixels. Radial ROI can be used to 
accurately measure MTF on PMMA phantom which has circular edge as long as the spatial resolution of the 
CT system does not pass the Nyquist frequency. 
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